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Abstract. Algorithm, based on explicit representations for analytic continuation of 
1^^ ' the T-matrix Faddeev components on unphysical sheets, is worked out for calculations 

0\ ' of resonances in the three-body quantum problem. According to the representations, 

CTs ■ poles of the T-matrix, scattering matrix and resolvent on unphysical sheets, inter- 

preted as resonances, coincide with those complex energy values where appropriate 
truncations of the scattering matrix have zero as eigenvalue. Scattering amplitudes 
on the physical sheet, necessary to construct scattering matrix, are calculated on the 
basis of the Faddeev differential equations. Effectiveness of the algorithm developed 
is demonstrated for example of searching for resonances in the system nnp and in a 
model three-boson system. 
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We make a numerical test of the approach proposed in Ref. ||l| to treat the three-body 
f~^ I resonances in the case of pairwise interactions falhng off in coordinate space not slower 

cn ■ than exponentially. This approach is based on a construction of explicit representations 

for the Faddeev components Map{z), a,f3 = 1,2,3, of the three-body T-matrix T{z) as 
r-^ ' well as for the scattering matrices S{z) and resolvent R{z) continued on unphysical sheets 

t^^ . of the energy z plane. For the sheets, the notation 11; is used with /, an enumerating multi- 

^1* I index (see |]l[| ) . The representations constructed demonstrate a structure of kernels of the 

above operators after continuation and give new capacities for analytical and numerical 
studies of three-body resonances. In particular the representations for analytic contin- 
ZJ ' nation M(z)\yi of the matrix M{z) = {Ma/3{z)} on the sheet Hi with details omitted 

^ . read 

M|n, = M|n„ - QJJ^ASf'jQi 
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The operator Qi{z) and the "transposed" one, QJiz) are obviously constructed of the 
matrix M{z) taken on the physical sheet Ho. A{z) is a number matrix, an entire function 
P^ \ of z G C. By Si{z) we understand a truncation (depending essentially on I) of the total 

three-body scattering matrix S{z). Operators J{z) and J'^{z) realize a restriction of 
kernels of the operators Qi{z) and Qi{z) on the energy shells, respectively, in first and 
last momentum variables. So that the products QJJ^ and JQi have half-on-shell kernels. 
Representations ||I| for analytical continuation of the three-body scattering matrices and 
resolvent follow immediately from the representations above for M{z)\yi ■ 

As follows from the representations constructed, the nontrivial (i.e. differing from the 
poles at the discrete spectrum eigenvalues of the three-body Hamiltonian) singularities of 
the T-matrix, scattering matrices and resolvent situated on an unphysical sheet 11;, are 
singularities of the inverse truncated scattering matrix Sf (z). Therefore, the resonances 
on the sheet 11; considered as poles of the T-matrix, scattering matrix and resolvent 
continued on 11; are those values of the energy z for which the matrix Si{z) has zero 

^LANL E-print pucl-th/9702037| . To appear in Proc. of Intern. Conf. on Computational Modelling 
and Computing in Physics (Dubna, September 16-21, 1996). 



Figure 1: Surface of the function |S'oi(z)| in the model system of three bosons with the 
nucleon masses. The potential V'^{r) is used with the barrier Vf, = 1.5 MeV. Position of 
the resonance Zres(3B) corresponds to the minimal (zero) value of |S'oi(z)|. 

as eigenvalue. Thereby, to search for the resonances situated on a certain unphysical 
sheet 11/, one can apply any method allowing to compute analytical continuation on the 
physical sheet of the elastic scattering, rearrangement or breakup amplitudes necessary 
for construction of respective Si{z). In particular such is the algorithm developed for 
(2 —^ 2, 3) processes in Ref. |Q] on the base of Faddeev differential formulation of the 
scattering problem in configuration space (see also |^ Q and Refs. therein). It is only 
necessary to go out in this formulation on the complex plane of z including the asymptotical 
boundary conditions. 

In the present work we utilize a code based on the algorithm H-H], for computations 
of s-state nnp resonances situated on the unphysical sheet 11(0,1) connected with physical 
one by crossing the spectral interval {E^, 0) between the deuteron energy z = E^ and 
breakup threshold z = 0. We solve the two-dimensional Faddeev integro-differential equa- 
tions Q with the (2 — > 2, 3) asymptotical boundary conditions |@]-0] at complex energies 
z and extract the truncated s-state scattering matrix Soi{z) = 1 -|- 2iao{z) with ao{z), 
the amplitude of elastic nd scattering continued on the physical sheet. When making a 
finite-difference approximation of the equations above in polar coordinates we take up to 
180 points of grids in both radial and angular variables, the cut-off radius being up to 
39 fm. As a A^A'-interaction, the Malfliet-Tjon potential MT I-III 0] is chosen. 

Firstly, we have checked a validity of the code finding the ^H bound-state energy Et 
as a pole of the function 501(2;). More precisely, the location of the Soi{z) pole was found 
as a root of the inverse amplitude l/aQ{z). Beginning from the grid dimension 80x80, we 
have obtained Et = —8.55 MeV. Hereafter all the energies are given with respect to the 



Figure 2: Trajectory of the resonance Zres(3B) on the sheet H/o^i) in the model system of 
three bosons with the nucleon masses. The potential V [r] is used. Values of the barrier 
Vb in MeV are given near the points marked on the curve. 

breakup threshold. Note that the value stated is in a good agreement with known results 
on Et in the MT HII model (see Ref. g]). 

Concerning the nnp resonances on the sheet H/o^i), we have inspected a domain of 
a range about 10 MeV in vicinity of segment [E^, 0] in the complex z plane. Especially 
carefully we studied a vicinity of the points 2: = — 1.5=b0.3 + i(0.6±0.3) MeV interpreted in 
Refs. 0, P as a location of an exited state energy of ^H. Unfortunately we have succeeded 
to find only one root z^es of the function Soi{z), corresponding to the known virtual state of 
the nnp system at total spin S=l/2. On a 180x180 grid, we have found Zres = —2.728 MeV 
i.e. it is situated 0.504 MeV to the left from the nd threshold Ed = -2.224 MeV (in the 
MT I-III model) . Note that the shift E^, — ^^res found from experimental data on nd 
scattering, is 0.515 MeV (see Ref. @)- Its value [^ computed in a separabilized MT I-III 
model on the base of the momentum space Faddeev equations, is equal to 0.502 MeV. As 
one could expect (see the data on three-nucleon resonances in ^), we have failed to find 
any resonances in the quartet state at L = as well as at L = 1. 

Also, we have studied a behavior of a resonance situated on the unphysical sheet 11(0,1) 
in a model three-body system including bosons with masses of the nucleon. As a pairwise 
interaction between the bosons we have used the Gauss-type potential supplied with an 
additional Gauss repulsive barrier term, 

V^{r) = Vo exp[-/xor2] + Vb exp[-/xfe(r - nf] 

where the values Vq = —55 MeV, fiQ = 0.2 fm~^, fib = 0.01 fm~^, r^, = 5 fm have been 
fixed while the barrier amplitude Vb varied. A resonance (with non-zero imaginary part) 
on the sheet 11(0,1) arises in the system concerned just due to the presence of the barrier 
term. Example of a surface of the |S'oi(2;)| function for the barrier amplitude V^ = 1.5 MeV 
is shown in Fig. || (for a 80x80 grid). A trajectory of the resonance 2:res(3B) (a zero of 
the function 5*01(2:)) is shown for the changing barrier Vb in Fig. ^. This trajectory was 
watched for the barrier Vb decreasing in the interval between 1.5 MeV and 0.85 MeV. When 
drawing the trajectory, we have used a 160 x 160 grid. It can be seen from Fig. |2| that the 
behavior of the resonance Zres(3B) is rather expected: with monotonously decreasing real 



Figure 3: Dependence of the "deuteron" energy E^ (curve 1) and real part of the resonance 
Zres(3B) (curve 2) on the barrier value V^. 

part, the imaginary part of the resonance changes also monotonously. For Vh > 0.95 MeV 
the energy Zi.es(3B) becomes real, turning into a discrete spectrum eigenvalue. In Fig. ^ 
we plot both the trajectories of the resonance real part Rezres(3B) and two-boson binding 
energy Ed- As one can see from this figure, the value of | Re2:res(3B)| increases more 
quickly than jE^^I, coinciding with \Ed\ at Vf, = 0.97 MeV. 

Trajectory of the resonance concerned in the lower complex half-plane is symmetric to 
the curve shown in Fig. |2| with respect to the real axis. Respective points, symmetric to 
those marked in Fig. |2|, correspond to the same values of V^. 
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